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New observations of ULX supershells, and their implications 
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New optical narrowband imaging observations of the fields of several ULXs are presented. Known supershell nebulae are 
associated with a number of these ULXs, which we detect in emission lines filters such as [S II], He II, [O II] and [O III]. 
New nebulae are discovered, which are candidate ULX-powered supershells. The morphologies and emission line fluxes 
of these nebulae could then be used to infer the properties of the emitting gas, which gives clues to the energizing source 
(photoionization and/or shock-excitation, both possibly from the ULX). Studies of supershells powered by ULXs can help 
to constrain the nature of ULXs themselves, such as the isotropy of the X-ray emission and the strength of their outflows. 
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1 Introduction 

Ultraluminous X-ray sources (ULXs) are objects emitting at 
X-ray luminosities exceeding ^ lO''^ erg s^^ (assuming the 
radiation is isotropic) that usually reside within galaxies but 
not at their centres. They are not considered to be supermas- 
sive black holes (e.g. active galactic nuclei; AGN), nor can 
they be stellar-mass compact objects (black holes or neu- 
tron stars) accreting at sub-Eddington rates from a compan- 
ion star (X-ray binaries; XBs). Their nature remains gener- 
ally elusive, although most are likely to be XBs harbouring 
black holes of tens (but not hundreds) of solar masses (Jin, 
Feng & Kaaret 2010) accreting in specific states which can 
appear super-Eddington (Gladstone, Roberts & Done 2009). 
Some of the most luminous ULXs (termed hyperluminous 
X-ray sources) may contain intermediate-mass black holes 
(hundreds of solar masses; e.g. Farrell et al. 2009). 

Unlike most XBs, but like many AGN, ULXs are often 
found to be associated with large, resolved structures. These 
are ULX nebulae (see e.g. PakuU & Grise 2008), ionized 
gas surrounding the ULX, typically found in optical line 
emission (many ULXs are also in proximity to star-forming 
regions; Swartz, Tennant & Soria 2009). Since ULXs are 
by definition strong X-ray/UV emitters, their radiation pho- 
toionizes the local interstellar medium (ISM), producing ob- 
servable nebulae tens of parsecs in size (e.g. the He III re- 
gion surrounding the ULX in Holmberg II; Kaaret, Ward, & 
Zezas 2004). Such nebulae can be used to infer the isotropic 
luminosity of the ULX, and hence test for relativistic beam- 
ing of its X-rays. The larger scale analogies are AGN, which 
heat and ionize the intergalactic medium, suppressing cool- 
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ing flows and galaxy formation. On smaller scales, an X-ray 
ionized nebula of radius ^ 5 pc is found associated with the 
black hole XB LMC X-1 (e.g. Cooke et al. 2008). 

Evidence is mounting that ULXs also release energy ki- 
netically via powerful, collimated jets/outflows, as is the 
case for both XBs and AGN. While the jets themselves are 
unlikely to be detected because their expected radio fluxes 
are fainter than current radio telescope detection limits 
(Kording, Colbert & Falcke 2005), ULX jets and/or winds 
have been inferred energetically via 'superbubbles' hundreds 
of parsecs in diameter and some bipolar in morphology, 
which are powered collisionally (e.g. Ramsey et al. 2006; 
Abolmasov et al. 2007; Abolmasov 2010). Historically, ra- 
dio lobes in FRJI radio galaxies are used as calorimeters to 
infer the kinetic power contained in the AGN jets fueling 
them. In the last few years these techniques have been suc- 
cessfully applied to the few known examples of XB jet-ISM 
interactions (e.g. Gallo et al. 2005; Pakull, Soria & Motch 
2010). Quite surprisingly, the results so far suggest that a 
large fraction of the accretion energy is channeled into the 
mechanical power in the jets (Dubner et al. 1998; Russell 
et al. 2007; Cooke et al. 2007; Sell et al. 2010; Soria et al. 
2010); a result that is a large step forward in our understand- 
ing of accretion physics in strong gravitational fields. The 
jet powers cannot be inferred from observations of the jets 
themselves since an uncertain fraction of the energy in the 
jets is radiated away (perhaps ^ 5%; Fender 2001). 

Unlike AGN, few large-scale jet-powered lobes/bubbles 
have been discovered so far associated with XBs (e.g. Mirabel 
et al. 1992; Miller-Jones et al. 2008), and indeed more can- 
didates may exist for ULXs. For XB jets, a high local mass 
density is required for a bow shock to form (Heinz 2002), 
the surface brightness of the radio lobes are typically be- 
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Table 1 Filters used. A is the central wavelength. The red- 
shifted Ha filter is denoted by Haz. 



Filter 




WFC 




DFOSC 




A (A) 


FWHM (A) 


A (A) 


FWHM (A) 


[O 11] 


3727 


100 






He 11 


4686 


100 






H/3 


4861 


170 






[0 111] 


5008 


100 


5010 


57 


Harris V 


5425 


975 






Hq 


6568 


95 


6562 


62 


Hqz 


6657 


79 






[S 11] 


6725 


80 


6727 


63 



low current detection limits (Kaiser et al. 2004), and a low 
space velocity of the system is needed so that the power is 
not dissipated over too large a volume (although such struc- 
tures could still be detected in some cases; Heinz et al. 2008; 
Wiersema et al. 2009). If the ratio between jet power and X- 
ray luminosity could be similar for ULXs and XBs, ULXs 
may be energetically more favourable for detections of jet- 
ISM interactions. In addition, their large extent of hundreds 
of parsecs would be hard to detect in our own Galaxy be- 
cause the nebulae would be too diffuse, spanning several 
degrees and requiring very wide field imaging. 

Observations of ULX nebulae can therefore provide a 
powerful tool to help uncover the nature of ULXs; placing 
constraints on their isotropic X-ray/UV luminosity and me- 
chanical output. Here we present optical narrowband obser- 
vations of the fields of several ULXs and identify a number 
of candidate ULX-powered nebulae. 

2 Observations 

Imaging of six galaxies containing a total of 17 ULXs (Liu 
& Mirabel 2005) were performed using two telescopes. We 
observed on three consecutive nights in April 2009 using 
the Wide Field Camera (WFC) on the 2.5-m Isaac Newton 
Telescope (INT) located at La Palma, Spain. The WFC has 
a field of view (FOV) of 34x34 arcmin. One target, NGC 
7793, was observed in September 2006 with the Danish 
1.5-m telescope at La Silla, Chile, equipped with the Dan- 
ish Faint Object Spectrograph and Camera (DFOSC; with 
a FOV of 13.7 x 13.7 arcmin). Some ULXs have known 
nebulae and some do not; for those which do we used many 
narrowband filters because the fluxes and emission line ra- 
tios can be used to infer the properties of the emitting gas. 
For instance, the [S II] doublet is most luminous for shock- 
ionized gas and the He II line is strongest for photoionized 
gas. For the WFC observations a redshifted narrowband Ha 
filter was used for continuum subtraction for filters centred 
in the red (the same exposure times were used); for filters 
centred at blue wavelengths the broadband Harris V filter 
was used (with a shorter exposure time). By subtracting the 
continuum emission, most of the galaxy light is removed 
and just the line emission is revealed. The filters used are 
listed in Table 1 . 




Fig. 1 Continuum-subtracted images of the region around 
NGC 5204 X-1 (north is up, east to the left), (a): Ha 
(Canada France Hawaii Telescope image; R-band is used 
for continuum subtraction; data first presented in Pakull & 
Mirioni 2002). (b-d): WFC images in [O III], He II and [S 
II], respectively. The FOV of each panel is 48 x 26 arcsec 
(- 1120 X 610 pc at a distance of 4.8 Mpc). The Chandra 
error circle of the position of the ULX (Goad et al. 2002) is 
shown in red. The inset in (c) is a radio image of the W50 
nebula distorted by flie jets of SS433 (Dubner et al. 1998), 
at the same size scale (assuming a distance of 5.5 kpc to 
SS433 and W50; Lockman, BlundeU & Goss 2007) show- 
ing that the NGC 5204 X-1 nebula is 2-3 times larger in 
length than W50. 

In addition, spectra of the optical counterpart of NGC 
5204 X-1 were obtained with the 8.1-m Gemini North tele- 
scope using the Gemini Multi-Object Spectrograph (GMOS) 
in July 2008. The slit used was 0.75 arcsec in width and 330 
arcsec in length, orientated north-south. The log of all ob- 
servations is presented in Table 2. 

Data reduction of the imaging data was performed us- 
ing the pipeline package THELl (Erben et al. 2005) for the 
WFC data and IRAF for the DFOSC data. The science im- 
ages were de-biased and flat-fielded using master bias and 
master flat frames. Positional calibration of the WFC data 
was achieved by matching several hundred stars in each ex- 
posure with those in the online MAST Guide Star Catalog 
(release 2.3) using THELl. The DFOSC images were aligned 
and stacked in IRAF to produce one deep image per filter. 
The GMOS spectra were reduced (de-biased, flat-fielded, 
sky subtracted, combined, flux calibrated and wavelength 
calibrated) using standard Gemini packages within IRAF . 

3 Results 

X-ray ionized nebulae are at most tens of parsecs from the 
ULXs energizing them and hence may be unresolved in our 
data (one example being NGC 5408 X-1; Kaaret & Corbel 
2009). However, ULX superbubbles that are energized colli- 
sionally are typically much larger, easily resolvable at these 
distances and may even appear detached from the ULXs 
since the evolution of jets causes underdense cavities sur- 
rounding XBs to develop (e.g. Hao & Zhang 2009). We do 
indeed find nebulae in proximity to several ULXs, some al- 
ready reported and some new. Here we present some inter- 
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Table 2 Log of observations. For filter wavelengths see Table 1. In the fourth column, 4x 10 refers to 4 exposures of 10 
minutes each. 



Date (UT) 


Instrument 


Galaxy 


Filters (exposure times in minutes) 


Airmass 


2006 September 1 


DFOSC 


NGC 7793 


Ha (4x 10), [O III] (3x 10), [S II] (3x10) 


1.00-1.18 


2008 July 3 


CMOS 


NGC 5204 


Spectra: wavelength range 3800-6000 A (4x 12) 


- 1.4 


2009 April 9 


WFC 


NGC 4861 


Ha (20+2), [S II] (20), Haz (20) 


1.01-1.04 


2009 April 9 


WFC 


NGC 4559 


Ha (20), [S II] (20), Haz (20) 


1.14-1.27 


2009 April 10 


WFC 


NGC 5204 


Ha (10), Haz (10), [S II] (10), V (2.5), [0 II] (10), H^ (10) 


1.20-1.30 


2009 April 10 


WFC 


NGC 4395 


Ha (20), Haz (20), [S II] (20) 


1.01-1.04 


2009 April 11 


WFC 


NGC 4559 


[O III] (20), [0 II] (20), He II (20) 


1.00-1.03 


2009 April 11 


WFC 


NGC 5204 


[S n] (10), [0 ni] (10), He II (lO), Haz (10) 


1.18-1.22 


2009 April 11 


WFC 


NGC 5194 


[S II] (10), Haz (10) 


1.16-1.18 



esting prehminary results from this study. In future work, 
we will use the emission line fluxes of the nebulae to infer 
their nature, estimate the photon and/or mechanical energy 
input and assess whether they could be powered by ULXs. 

NGC 5204 X-1 is a ULX with a variable X-ray lumi- 
nosity of 1-6 xlO^^ erg s~^ (0.5-8 keV) residing in a 
warped spiral galaxy at a distance of 4.8 Mpc (e.g. Roberts 
et al. 2006). Roberts et al. (2001) detected ionized gas within 
100 pc of the ULX in Ha, H/3, [S II] and [O III], with an 
apparent cavity on one side; Abolmasov et al. (2007) found 
several emission lines in the optical spectra; and a 360 pc 
diameter circular 'bubble' with an even larger, apparently 
bipolar structure was seen in an Ha image by Pakull & Miri- 
oni (2002). In Fig. 1 we present some of our continuum sub- 
tracted WFC images of the environment of this ULX (in Fig. 
la the high resolution Ha image first presented in Pakull 
& Mirioni 2002 is shown). Regions of the circular bubble 
are detected in [O III] (Fig. lb), suggesting that this may 
be a shock wave travelling at > 100 km s^^ (e.g. Allen 
et al. 2008). Abolmasov et al. (2007) also showed that a 
nebula close to the ULX could be energized collisionally 
and/or photoionized, from spectra. Since this nebula is part 
of the circular bubble, the morphology suggests a shock 
wave (which could also be partly photoionized). The two 
ionized regions to the north and south of the bubble are de- 
tected in [O 111] (Fig. lb), [S 11] (Fig. Id), [O 11] and H/3 
(not shown) and faintly in He II (Fig. Ic). 

Numerous emission lines from the nebula are visible in 
the GJVIOS spectrum of NGC 5204 X-1 (Gladstone 2009). 
The sUt is orientated north-south, and emission from the 
ULX optical counterpart, the circular bubble and large scale 
bipolar structure are all detected. By comparing the de- red- 
dened [O II] / H/3 and [Ne III] / HP emission Une ratios 
with models for radiative shock waves (Allen et al. 2008), 
initial results suggest that much of the gas in the large-scale 
nebula is travelling at velocities of the order ~ 200 km s~^ . 
The non-detection of He 11 in the GMOS spectrum and only 
faint detection in some regions in the WFC image imply that 
photoionization does not play a large role in the energetics 
in most regions of the superbubble; the bright [O 11] and 
[S 11] seen in the WFC images of the bipolar structure im- 
ply it could be coUisionally energized. The brightest parts 
of the whole structure appear to form a straight line that al- 



most passes through the position of the ULX (most visible 
here in the [O III] image). These preliminary results suggest 
this system could represent the first known ULX analogue 
of the large-scale SS433 / W50 complex (e.g. Dubner et al. 
1998), in which the jets of SS433 are distorting and energiz- 
ing the supernova remnant W50. W50 is 2-3 times smaller 
in length than the NGC 5204 X-1 superbubble (see inset in 
Fig. Ic); SS433 also has a lower X-ray luminosity. 

NGC 4559 is a late-type spiral galaxy at a distance of 
~ 9.7 IVIpc (Sanders et al. 2003) that contains two ULXs; 
XIO (Lx ~ 1 X 10*0 erg s'^ 0.3-10 keV) residing close 
to the centre of the galaxy (but not consistent with its nu- 
cleus) and X7 (Lx ~ 2-3 x 10*" erg s-^; 0.3-10 keV) in 
an isolated, outer region of the galaxy (Cropper et al. 2004). 
X7 is located within a large star-forming complex, studied 
extensively by Soria et al. (2005). The authors show that a 
ring-shaped Ha structure is coincident with the star-forming 
regions and cannot be ionized (solely) by X-rays from the 
ULX. In Fig. 2 the fields of these two ULXs are shown 
in continuum-subtracted emission Une images taken with 
the WFC. We find that the Ha ring surrounding X7 is also 
bright in [0 11] and [O III] (upper panels). Fainter detections 
are also made in [S 11] and He 11 (not shown here). The line 
emission likely originates in gas ionized by the population 
of young massive stars, but we can use the emission line 
ratios to assess the possibiUty of photoionization or shock 
excitation from the ULX (in a follow-up work). An emis- 
sion line nebula hes a few hundred parsecs north-east of 
the position of XIO (Fig. 2, lower panels). The western re- 
gion of the nebula, closer to the ULX is brighter in [O 11], 
[O III] and He II, whereas the eastern region is brighter in 
[S II]. This is a new candidate ULX superbubble, but fur- 
ther analysis is required to test whether this nebula could be 
energized by the ULX. 

4 Discussion and conclusions 

We have introduced a study into the ionized environments 
of a selection of ULXs, and give some preliminary results. 
Narrowband images were obtained using two telescopes in 
many emission line filters sensitive to photoionized and shock- 
ionized gas. Candidate ULX-powered supersheU nebulae are 
identified, including known and new candidates. In a future 
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Fig. 2 Continuum-subtracted images of the region around 
NGC 4559 X7 (upper panels) and XIO (lower panels). Each 
panel is 30 x 30 arcsec (« 1.4 x 1.4 kpc at a distance of 
9.7 Mpc); north is up, east to the left. The red circles indi- 
cate the X-ray positions of the two ULXs (the error circles 
are smaller than the ones shown). We measure the positions 
from archival Chandra data using the source detection pro- 
gram celldetect (Calderwood et al. 2001); these agree with 
the positions published in Cropper et al. (2004). 

work we will focus on individual ULXs and use the mor- 
phologies and line ratios to infer the likely energy input to 
these structures from the ULXs. 

Although the environment of ULXs has only began to 
be studied in detail in the last decade, the number of re- 
gions in the environment likely energized, ionized or dis- 
rupted by ULXs has already exceeded the number of known 
cases for canonical XBs. This may be due to a combination 
of angular size scale issues (such diffuse structures in our 
galaxy are harder to identify since they would extend ar- 
cminutes or degrees across), lower extinction for extragalac- 
tic sources, the more luminous nature of ULXs compared 
to XBs and likely denser regions of the ISM around ULXs 
(which have young, high-mass companion stars) compared 
to low-mass XBs (which have older, low-mass stars). ULX 
nebulae may be abundant, and observable at many wave- 
lengths with current facilities. Systematic, unbiased sam- 
ples may prove pivotal in advancing our understanding of 
ULXs and the feedback they provide in heating the ISM 
in galaxies. The largest superbubbles appear to be at least 
partly jet-powered. For XBs, the kinetic energy of jets is 
hard to measure, but is a remarkably large fraction of the ac- 
cretion energy in the few cases where estimates have been 
possible. ULXs may therefore provide the best opportuni- 
ties to measure the energy contained in the jets produced 
by stellar-mass (and possibly even intermediate-mass) black 
holes. 
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